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Abstract 

We discuss recent results in the study of the evolution of strongly coupled field 
theories in the presence of time dependent couplings using the holographic corre- 
spondence. The aim is to understand (i) thermalization and (ii) universal behavior 
when the coupling crosses a critical point. Our emphasis is on situations where a 
subset of bulk fields can be treated in a probe approximation. We consider two 
different setups. In the first, defect conformal field theories are described by probe 
branes in AdS space-times, and an initial vacuum state evolves due to a time de- 
pendent coupling in the probe sector. While a black hole formation is invisible in 
this approximation, we show that thermalization can nevertheless happen - this is 
signalled by formation of apparent horizons on the brane worldvolume. In the sec- 
ond setup, we consider a probe bulk scalar field in the background of a AdS black 
brane. In equilibrium, this system undergoes a critical phase transition at some 
temperature when the source for the dual operator vanishes. For a time dependent 
source which goes across this critical point, we show that a zero mode of the bulk 
field dominates the dynamics and leads to scaling behavior of the order parameter 
as a function of the rate of change. 
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1 Introduction 



The problem of quantum quench, i.e. the response of a system to a time dependent 
coupling, has recently attracted a lot of attention in several areas of many-body physics, 
particularly because of progress in cold atom experiments [I], [2], [3]. This problem is 
interesting for at least two reasons. The first relates to the question of thermalization. 
Suppose we start with the ground state and then turn on a time dependent coupling 
which again approaches a constant at late times. Does the system evolve into some kind 
of steady state ? If so, is the state "thermal" in any sense ? 

The second question deals with the situation where the quench takes place across a 
value of the parameter where there is an equilibriium critical point. In this case there 
is some evidence that the time evolution carries some universal features of the critical 
point. Assuming that there is only one scale which governs the behavior of the system in 
the critical region, one can, e.g. derive scaling properties of one point functions following 
adaptations of the early work of Kibble and Zurek [HE]. Suppose the coupling approaches 
the critical coupling linearly, i.e. 

{g-g c )~vt (1.1) 

These arguments then show that the one point function of an operator with conformal 
dimension x at the critical point has a scaling behavior [1], 

< 0{t) >~ (v)^TiF(tv^i) (1.2) 

where z is the dynamical critical exponent and v is the correlation length exponent. 
Another manifestation of such universal behavior appears in 1 + 1 dimensional field the- 
ories which are quenched suddenly to a critical point. In this case, powerful methods of 
boundary conformal field theory can be used to obtain the time dependence of correlation 
functions [2j|3] . For example the one point function of a generic operator with conformal 
dimension x behaves as 

<£>(t)>~e"^ (1.3) 

where To denotes a length scale which characterizes the initial theory away from criticality. 
tq is of course not universal, so neither is the relaxation time r = 2to/(itx). However, 
the ratio of the relaxation times for two different operators 0\ and O2 with conformal 
dimensions X\ and x^ is universal, 

^ = ^ (1.4) 

r 2 Xi 

Unlike equilibrium critical phenomena there is no general theoretical framework to under- 
stand such scaling relations. In particular, there are very few theoretical tools available 
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to study such systems when they are strongly coupled. It is, therefore, natural to explore 
if the AdS/CFT correspondence [6] - [9] is useful in this problem. 

In this contribution we will summarize results obtained in this approach. The work on 
thermalization is in collaboration with Tatsuma Nishioka and Tadashi Takayanagi |10j . 
The work on quench across critical points is in in collaboration with Pallab Basu 

In the AdS/CFT correspondence, couplings of the boundary field theory are boundary 
values of a bulk field. In the regime where supergravity is valid, the problem of quantum 
quench then reduces to a classical problem with given initial and boundary conditions. 
This problem has been studied when the time dependent coupling is the boundary metric 
or the gauge theory coupling (i.e. the boundary value of the dilaton). Suppose this 
coupling is a constant in the far past and future, and has a smooth time dependent profile 
at intermediate times. In the bulk description this corresponds to a disturbance created 
on the boundary which propagates in the bulk. Under suitable conditions, this leads to 
black hole formation in the bulk [T2T[Ti] . The correlators at future time would then be 
thermal with a temperature characterized by the Hawking temperature. The time scale 
after which this happens depends on the nature of the correlators, but turns out to be 
always smaller than what one would expect from a conformally invariant system evolving 
to a thermal state. Thus, in this case thermalization of the field theory is signalled by 
black hole formation. For homogeneous planar collapse (i.e. a space independent coupling 
in Poincare patch), a black hole is always formed. In the case of homogeneous collapse in 
global AdS, a black hole is formed when the rate of change is fast enough compared to 
the scale set by the radius of the sphere on which the boundary theory lives. 

In other situations, e.g. a slow variation of the coupling for global AdS, a black hole 
is not formed so long as the supergravity approximation is valid. Rather, if the coupling 
becomes weak at some time, the bulk string frame curvature grows large, leading to a 
breakdown of the supergravity approximation - thus mimicking a space-like singularity 
[15]. For the case of a slow variation of the coupling it turns out that the gauge theory 
remains well defined and may be used to show that a smooth passage through this region 
of small coupling is possible without formation of a large black hole. Related scenarios 
appear in [16] and [17] . 

One of our main aims is to study quantum quench across critical points. Many such 
critical points can be studied in setups where a subset of the bulk fields can be treated 
in a probe approximation. In this approximation, these probe fields provide the essential 
physics and their backreaction to the background gravity can be ignored, typically sup- 
pressed by 1/N. This motivates us to study the problem of quantum quench in situations 
where such a probe approximation is valid. 
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We will first consider defect field theories which arise as dual descriptions of a set 
of probe branes in the AdS x S bulk [18J. This approach has been used extensively to 
study flavor physics, as well as models with possible applications to condensed matter 
systems. The nice feature of this approach is that the boundary field theory is known, 
though they typically have supersymmetry. In this case a quantum quench of couplings 
in this subsector becomes a classical motion of these probe branes, with specified time 
dependent boundary conditions at the AdS boundary [10]. We investigate the question 
of thermalization in this context. Since the background geometry is unchanged in this 
approximation any black hole which is formed due to the quench is not visible. We will 
find that thermalization is nevetheless visible - this manifests itself as the formation of an 
apparent horizon on the brane worldvolume. 

In the second setup we consider a "bottom- up" bulk theory of gravity with a neutral 
scalar field [TJJ], where one writes down a bulk theory and assume that there is some 
dual field theory on the boundary. In this specific instance, the background is a AdS± 
charged black brane and the coupling of the scalar is large, so that its backreaction 
to the geometry is small. When the mass of the scalar lies in the range — § < m 2 < 
— | there is always a critical temperature below which the trivial solution is unstable. 
Equivalently, for a given temperature, whenever the mass is below a certain value, the 
trivial solution is unstable. In this regime there is a new nontrivial static, stable solution 
whose "non-normalizable" part vanishes. This means that in the dual theory there is a 
new phase where the expectation value of the operator dual to this scalar is non-zero, even 
in the absence of any source. For any nonzero temperature there is a continuous phase 
transition with mean field exponents at the critical mass. At zero temperature (i.e. when 
the background is an extremal brane) the transition persists, but is of the Berezinskii- 
Kosterlitz-Thouless type. This setup is similar to that of holographic superconductors 
[20H22] and has been proposed as models for antiferromagentic transitions. 

We consider quench across this critical point by working at the critical mass, but 
turning on a time dependent source which crosses zero (i.e. the critical point) at some 
time [H] . We show that the dynamics of the bulk scalar is dominated by a zero mode of 
the radial operator in the critical region when the rate of change of the source is small, 
This leads to a Landau-Ginsburg type dynamics with dynamical critical exponent z = 2, 
and a resulting scaling behavior of the order parameter. 
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2 Probe Branes and Thermalization 



Probe branes in the bulk of AdS have been used to introduce flavor in the standard 
AdS/CFT correspondence. Consider for concreteness AdS§ x S 5 whose dual is M = 4 
super- Yang-Mills in 3 + 1 dimensions with gauge group SU(N C ). Let us introduce Nf Dp 
branes which wrap a AdS m x S p+1 ~ m . Possible supersymmetric wrappings are summarized 
in Table (0. 



Table 1: Probe Branes in AdS$ x S 5 . 



Brane 


Wrapping 


Dual Theory 


Dl 


AdS 2 


+ 1 dim 


D3 


AdS 3 x S 1 


1 + 1 dim 


D5 


AdS A x S 2 


2 + 1 dim 


D7 


AdS 5 x S 3 


3 + 1 dim 



The Dp branes give rise to new hypermultiplet fields. These live on the intersection of 
the Dp branes with the D3 branes which gave rise to the AdS§ x S 5 geometry. From the 
point of view of the N = 4 theory the hypermultiplets live on a lower dimensional defect. 
In the strong coupling regime, the bulk theory is the original supergravity together with 
the action of branes coupled to it. 

In the limit of Nf <C iV c the backreaction of the probe branes on the background 
AdS§ x S 5 geometry can be ignored and the entire bulk theory is given by the action of 
these branes moving in the fixed background geometry. We will take the brane action to 
be of the DBI type. In the dual theory this means we can consider the defect field theory 
by itself, and ignore the effect of hypermultiplet loops. 

As is standard in the AdS/CFT correspondence, the boundary values of the DBI fields 
are identified with sources for the dual operators in the dual field theory. Consider for 
example the case of a D5 brane. Let us write the AdS$ x S 5 metric in the form 

ds 2 = (y 2 + r 2 ) \-dt 2 + dx\ + dx\ + dxl] + -^-^ \dr 2 + r 2 dn\ + dy 2 + y 2 d(n' 2 ) 2 } (2.5) 

yl _|_ r Z 

The D5 brane is wrapped along £ a = (t, r, Q2, Xi, x^)- The fields in the DBI action 
are ?/(£), ^(Oj x 3 (£), which are the transverse coordinates to the brane. The value of 
y(r = 00) is then the mass of the hypermultiplet fields coming from (3, 5) open strings 
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joining the D5 brane with the stack of N c three branes which produce the background 
geometry. Thus a time dependent boundary value of y is a time dependent mass for the 
hypermultiplets. 

Therefore a quantum quench in this dual theory may be implemented simply by pro- 
viding a time dependent boundary condition for the DBI field. However the DBI fields 
are the transverse coordinates of the branes - so this corresponds to a motion of the edge 
of the brane. This disturbance sets up a wave along the brane and therefore correponds 
to an excited state of the defect field theory. Our aim is to figure out the nature of this 
state at late times. 

In the full theory, such a disturbance would lead to a deformation of the background 
geometry and possibly lead to black hole formation, which would appear as thermalization 
in the boundary field theory. We want to explore if any signature of thermalization remains 
in the probe approximation. The following sections summarize some salient points of work 
with Tatsuma Nishioka and Tadashi Takayanagi jlOj . 



2.1 Rotating Dl branes 



The essential physics is in fact apparent in the simplest example - Dl branes in AdS^ x S 5 . 
For this purpose it is convenient to write the AdS$ x S 5 metric as 

ds 2 = 2drdv - f{r)dv 2 + r 2 ds 2 + {d6 2 + sin 2 6d<p 2 + cos 2 6dQ 2 3 ) (2.6) 

If we are using the Poincare patch, f(r) = r 2 and ds 2 is the flat metric on R 3 , while 
in the global patch /(r) = 1 + r 2 and ds 2 is the round metric on S 3 . We have used 
Eddington-Finkelstein coordinates in the AdS§ part. The Dl brane is along (r, v) and 
its action is given by the standard DBI action obtained from the induced metric. The 
dynamical fields on the brane are d(r,v),(p(r,v),Qs(r,v) and the coordinates contained 
in ds 2 . It is clear from the symmetries that one can have a class of solutions of the form 

V(r,v), Q = \ (2.7) 

with all the other coordinates held constant. The equations of motion which follow from 
the DBI action is best written by using the advanced EF coordinate u = v — 2 J j^- as 
well as v 

o o 2 fd u (pd v (p\ 2 ( d u <fd v <f\ 

o u d v ip + -d v (pd u I I + -d u (pd v I I = , (2.8) 

where 

L = 1 - -^d u <pd v <p . (2.9) 

fin 
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It is easy to see that any cp which satisfies either d u ip = or d v (p = is a solution of (12. 8p . 
In particular, a solution which is a function of v alone represents the retarded effect of a 
boundary value of <p. 

The induced metric produced by such a retarded solution is given by 



ds 2 



ind 



-f(r)dudv + (d v <p) 2 dv 2 = Idrdv - [f(r) - {d v ip) 2 ]dv 7 



(2.10) 



This is a two-dimensional AdS Vaidya metric which has an apparent horizon at /(r) = 
(d v (p) 2 , provided this equation has a solution for real r. In the Poincare metric there is 
always a solution for real r, while in global coordinates, this is not guaranteed. 

Depending on the profile of (p(v), the apparent horizon may or may not develop into 
an event horizon. An example where it does is given by the profile 



(p(v) = (fo(v + — log cosh(A;t>)) 
k 



(2.11) 



which leads to the following equation for the location of the apparent horizon for the 
Poincare patch 

r = <p (l + tanh(kv)) (2.12) 

This asymptotes to an event horizon at r = 2<^o- This profile represents a Dl brane which 
starts from rest and spins with an increasing spin, asymptoting to a constant rotation 
rate. The function (12.111) and the location of the apparent horizon is shown in Figure ([1]) 
and (ED 




Figure 1: The profile (12. lip as a 
function of v 



Figure 2: Location of the appar- 
ent horizon, v as a function of 
r from equation (I2.12p . At late 
time this becomes an event hori- 
zon 



For strings which eventually stop spinning, the apparent horizon does not develop into 
an event horizon, but recedes back to r = 0. For example if 



ip(v) = <po(l + tanh(/o> )) 



(2.13) 
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the apparent horizon is located at 



cosh (kv) 



(2.14) 



The function (I2.13P and the location of the apparent horizon is shown in Figure ([3]) and 





Figure 3: The profile ( 12,131) 
function of v 



Figure 4: Location of the appar- 
ent horizon, v as a function of 
r from equation (12.141) . The ap- 
parent horizon now recedes back 
to r = 



In global AdS the equation which determines the apparent horizon is given by 1 + r 2 = 
(d v ip) 2 which does not always have a solution for real r. For example, a brane whose end 
point is uniformly rotating has ip(v) = u>v - this would lead to an event horizon only when 
oj > 1. 

The Poincare patch solutions represent injection of energy from the boundary, which 
flows into the Poincare horizon. In the global solutions, the energy flows from a point of 
the boundary to the antipodal point of the S 3 . 

Fluctuations of the brane around this classical solution will feel the effect of an ap- 
parent horizon on the worldsheet. Let us choose a static gauge where the worldsheet 
coordinates are identified with two of the space-time coordinates £ a = x a , a = 0, 1. The 
transverse coordinates are x 1 , 1 = 2 ■ • • 9. The AdS x S metric can be then written as 



ds 2 = g ab (x a , x I )dx a dx h + G u (x a , x^dx^x 3 
Expanding around a classical solution 

x \ x a ) = x I Q {x a )+y I {x a ) , 
This leads to the following action for quadratic fluctuations 



So 



^V = ^^G I j{^,x I )9ay I d b y J . 



(2.15) 



(2.16) 



(2.17) 
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where 7g 6 denotes the induced metric due to the background solution x^. In particular, the 
fluctuations of (p, 9 i.e. all fluctuations in S q directions are minimally coupled massless 
scalars on the worldsheet, while the fluctuations of the boundary gauge theory spatial 
directions x l ,i = 1- ■ - 3 have an additional factor of r 2 coming from the fact that = 
r 2 Sij along these directions. 

It is well known that fields which live on a space-time with an apparent horizon 
behave approximately thermally if the apparent horizon lasts long enough [23] . While the 
standard derivation of Hawking radiation assumes the presence of an event horizon, the 
essential physics is the large redshift, which is present near an apparent horizon as well. 
In our case, profiles like (12. lip lead to exact thermality at late times since the apparent 
horizon evolves into an event horizon. On the other hand profiles like (I2.13P lead to an 
effective "time dependent temperature" in the dual theory, which of course makes sense 
when the time variation is slow enough. 

The thermal nature of the state produced by time dependence becomes clear from a 
calculation of the fluctuation of the end-point of the string. In |24j it has been shown 
that the fluctuations of a string suspended from the horizon of a AdS black brane ended 
at a flavor D-brane near the boundary of AdS are dual to Brownian motion of the corre- 
sponding quark in the hot M = 4 gauge theory. In this case the bulk black brane metric 
induces a worldsheet metric which has a horizon. The fluctuations then reflect Hawking 
radiation from the worldsheet horizon. 

In the D-brane solutions considered above, the bulk metric has no horizon. However 
due to the motion of the D-brane, the induced metric on the worldvolume can develop a 
horizon. Since the fluctuations of [24] comes purely from properties of the induced metric 
it is natural to expect that a similar phenomenon appears in our case. 

The result of this calculation for fluctuations in the p direction is 



((Ay*(t-t>)T) ~ n(t-t')«(3 



<(Ay*(f-f)) 2 > ~ *^-^log[27r(f-0/)8], tt(* -*')»/?, (2.18) 

while for fluctuations in the x direction are 

([Ay\t - 1')] 2 ) ~ t«m/3 2 , 

mp 

~ /3\t-t'\ , t >m/3 2 . (2.19) 

The rotating D\ brane solution corresponds to a time dependent coupling in the 
iV = 4 theory coupled to hypermultiplets living on the zero dimensional defect. The Dl- 
D3 system is 1/4-BPS and the D1-D3 open strings lead to the two complex scalars (Q, Q) 



8 



of hyper mult iplets which belong to the fundamental and anti-fundamental representations 
of the color SU (N) gauge group. Let us express the three complex adjoint scalar fields in 
the M = 4 super Yang-Mills by ($1, $2, ^3)- These correspond to cartesian coordinates 
in the transverse C 3 composed of (r, Q 5 ) where f2 5 represents the 5-sphere. We choose 
$3 such that its phase rotation describes the one in the ip direction and that 9 = 7r/2 is 
equivalent to $1 = $ 2 = 0. The time dependent coupling term which corresponds to a 
uniformly rotation Dl-brane is given by 

J dt [Q [Im($ 3 e-^)] 2 Q + Q [Im($ 3 e^)] 2 0\ • (2-20) 

The justification for this is given in [10] . For non-uniform rotation with a profile <p(v) the 
exponential factors are simply replaced by e ±l</3( ^. 

Thus, from the boundary theory point of view we have a time dependent coupling - 
this leads to thermalization. Note that this is thermalization of only the hypermultiplet 
sector - the vector multiplet sector is unchanged in the lowest order of this approximation. 



2.2 Higher dimensional probes 



Thermalization in higher dimensional field theories can be also investigated by considering 
higher dimensional branes in the bulk of AdS x S. It is possible to construct uniformly 
rotating D7 and D5 branes using a combination of analytic and numerical methods. 
While these D5 and D7 solutions have been obtained earlier in [25], the implications to 
thermalization was not realized. Quench-like solutions (i.e. solutions where the rotation 
vanishes in the asymptotic past and the asymptotic future) can be obtained numerically 
as well [10J. 

Rotating D5 branes are dual to a time dependent mass of the hypermultiplets in the 
dual 2+1 dimensional defect field theory. Rotating D7 branes are dual to a time dependent 
phase of the mass of fermions in the hypermultiplet, as well as a time dependent bosonic 
potential. The essential physics is similar to the Dl brane descibed in the previous 
subsection, viz. an apparent horizon is formed and the fluctuations respond in a thermal 
fashion. 

An additional signature of dissipation appears when we consider probe D3 branes 
obtained by performing T-duality on the Dl brane solution described above along x± and 
X2 directions. In this case the worldvolume is a 3 + 1 dimensional theory and the dual 
defect field theory is 2 + 1 dimensions. It turns out that in this case one can turn on 
a background electric field on the worldvolume of a uniformly rotating D3 brane. This 
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means we are turning on a chemical potential and a charge density for the corresponding 
global charge in the dual field theory, in addition to a time dependent coupling. The 
fluctuations of the gauge field around this background now react to the apparent horizon 
of the induced metric. This leads to an electrical conductivity cr(u) , whose behavior as a 
function of the frequency v is quite similar to Drude theory at low frequencies. However 
the real part of o[y) approaches a constant at large frequency, as is typical in a 2+1 
dimensional critical theories. This is shown in Figures fl5]) and (J6]) 




Figure 5: Real part of o{y). Figure 6: Imaginary part of o{y). 



2.3 Time Dependent Chemical Potential 



An interesting example of the process of thermalization due to formation of an apparent 
horizon on the worldvolume concerns the thermalization of the meson sector of N = 4 
Yang-Mills theory due to a time dependent chemical potential [26J. As usual, quarks 
are introduced by placing a set of D7 branes and a chemical potential corresponds to a 
worldvolume electric field. This is made time dependent by coupling to a time dependent 
external current, i.e. by injecting quarks from outside. As in the previous examples this 
results in the formation of an apparent horizon with a characteristic temperature. For 
other effects of this phenomenon see [27] . 

Finally the formation of apparent horizons in these examples is similar to acceleration 
horizons on worldvolumes discussed in and related phenomena have been studied 
in m\. 
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3 Quench across a Holographic Critical Point 

As we discussed above, quantum quench is particularly interesting when the time depen- 
dent coupling crosses a critical point. The probe approximation has turned out to be 
quite useful for studying holographic critical points. Such phase transitions are known to 
occur for various probe branes [301432] . 

Another class of probe fields appear in discussions of holographic superconductors 
[2"0H2"2"] . This setup consists of a charged scalar field in the presence of a charged black 
brane. When the gauge coupling is large, the backreaction of the scalar and the gauge field 
to the background geometry can be ignored. In this case, for a given mass of the scalar 
field there is a critical temperature below which the scalar condenses - this is interpreted 
as superfluidity in the boundary theory. The phase transition between the ordered and 
disordered phases is a critical point. 

An even simpler setting consists of a neutral scalar field with a quartic coupling which is 
large enough to ensure that a probe approximation is reliable. As shown in [19] for suitable 
values of the parameters, this model displays a critical point of the type encountered in 
antiferromagentic phase transitions. In this section I will describe a recent study of 
quantum quench in this model in collaboration with Pallab Basu [11] . 

3.1 The equlibrium phase transition 

The model of [T9] has a neutral scalar field <p(t, r, x) in the background of a charged AdS^ 
black brane. The lagrangian is given by 




(3.21) 



The background metric is given by (in Racls = 1 units) 



ds 2 = [-r 2 f(r)dt 2 + r 2 dx 2 ] + 



dr 2 



(3.22) 



r 



,2 /(r) 



where 




(3.23) 



The associated Hawking temperature is then given by 




(3.24) 
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In the following we will replace r — > rr^. This means all dimensional quantities are 
expressed in units of r . 

In the limit of large A the field can be regarded as a probe field. In [19] it was shown 
that when the mass lies in the range 

4 2 v ; 

there is a critical phase transition at some value of T = T c (m) when the source to the 
dual operator vanishes. Conversely, for a given T there is a value of m 2 = m 2 where the 
theory is critical. 

The upper limit in (I3.25f) is the BF bound for the near-horizon AdS2 geometry which 
appears in the extremal (rj = 0) metric. (Note that the AdS scale for this infrared AdS2 
is given by 1/y/E in our units). The lower bound is the BF bound for the asymptotic 
AdS 4. Field configurations which are translationally invariant in the x directions satisfy 
the equations of motion 

4[-7TT 9 t + d r (r 4 f(r)d r )]4> - m 2 - 3 = (3.26) 
r z j[r) 

Near the ArfiSVboundary the asymptotic behavior of the solution to the linearized equation 
is of the form[j 

0(r) = J(t)r- A - [1 + 0(l/r 2 )] + < O > (t)r~ A + [1 + 0(l/r 2 )] (3.27) 
where A is given by 

A± = \ ± \J™ 2 + \ (3-28) 

In the range of masses of interest, both the solutions are normalizable, so that there is 
a choice of quantization. The standard quantization considers the coefficient J(t) as the 
source in the dual field theory and B(t) then gives the expectation value of the dual 
operator. In the alternative quantization the expectation and source change the role. 

Consider first the linearized problem, ignoring the cubic term. By a standard change 
of coordinates to tortoise coordinates p and a field redefinition to \i 

"P=-4r, «r,«) = *M (3.29) 

The horizon is then at p = 00 and the boundary is at p = 0. At the linearized level, the 
equation (I3.26P becomes 

- dh = -dlx + V (p)x = V pX (3.30) 



3 When we turn on J(t), it is a valid concern whether we will be able to neglect the non-linear term 
near the boundary. This can be done as long as A > or m 2 < 0. 
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with 

Vo(p) = r 2 f(r)[(m 2 + 2) - ^ + 1±^] (3.31) 

where in Vo(p) we need to express r in terms of p using (13.291) . 

For solutions of the type \ ~ e~ luJt , equation ( 13 .311) is a Schrodinger problem in a 
potential Vo(p). The potential goes to zero at the horizon p = oo and behaves as < " m p t 2 ^ 
near the boundary p = 0. Note that for a brane background at any finite temperature, 
f(r) ~ (r— 1) near the horizon, while p ~ — log(r — 1) so that Vo ~ e _p as we approach the 
horizon. A typical Vo(p) for m 2 > —2 is shown in Figure([7]). In contrast, for the extremal 



v_o 




Figure 7: The potential V (p) 

background f(r) ~ (r — l) 2 while p ~ l/(r — 1) so that Vq ~ 1/p 2 - This makes the 
analysis for the extremal background rather subtle. We will work with the non-extremal 
case. 

In [19] it was shown that when m 2 is below a critical value, m 2 there are bound states 
of this Schrodinger problem, showing that the trivial = is unstable. At m 2 = m 2 a zero 
energy bound state appears, which vanishes in an appropriate fashion at the boundary and 
is in addition regular at the horizon. In the complex frequency plane some quasinormal 
mode(s) hit the origin at m 2 = m 2 c . This critical mass is m 2 = — | when r\ — 1 and 
decreases with decreasing rj or increasing temperature. 

For m 2 < m 2 there is a stable nontrivial static solution 0o( r ) of the full nonlinear 
equation of motion with the condition that J = 0. This means that in the conventional 
quantization, the expectation value of the dual operator is nonzero even in the absence 
of a source, i.e. the field condenses. The the critical point is at m 2 = m 2 and J = 0. 
Similarly there is a nontrivial solution with B = which means that there is a condensate 
in the alternative quantization as well. 

The critical point has mean field exponents at any finite T. If the operator dual to 
the field <fi is O and the source is J then an analysis identical to that presented in [19] 
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leads to (for m 2 — m 2 — >■ + ) 

< O > J=0 ~ {m 2 -m 2 c f/ 2 ^^\j=o ~ (m'-ml)- 1 < O >„ ~ J 1/3 (3.32) 

Exactly at zero temperature the phase transition is of BKT type and the order parameter 
depends exponentially 

7Ta/6 



< O >j =0 ~ exp 

2 V m c ~ 

3.2 Quenching across the Critical Point 



m 2 



(3.33) 



In the following, we will study quench across this critical point by considering a time 
dependent source J(t) which asymptotes to constant values at early and late times and 
crosses zero at some time, e.g. 

J{t) = J tanh«) (3.34) 

In the conventional quantization this means that in the dual boundary field theory, we 
have a source coupling to the operator dual to 0. For any static J ^ we of course have 
a nontrivial <p(r) and hence a nonzero < O >. Our first aim is to get some insight into 
the time dependence of < 0(t) > when we have a nontrivial J(t). 



3.3 Breakdown of Adiabaticity 



If the time dependence is slow enough one would expect that far away from the critical 
point the dynamics is adiabatic, while near the critical point adiabaticity should break 
down. It is instructive to examine the way this happens. 

It is well known that to study low frequency modes in the background of a black brane 
it is convenient to use ingoing Edddington-Finkelstein coordinates, 

u = t — p, p (3.35) 

where p is defined in (13.291) . In terms of these coordinates the equation of motion [3.26j) 
becomes 

- 2d u d pX = -d 2 pX + V(j>, x). (3.36) 

where 

V(p,x) = V (p) X + f(r)x 3 - (3.37) 
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This equation has to be solved with the boundary condition that the field is regular at the 
horizon, which at the linearized level is equivalent to requiring that the waves are purely 
ingoing at the horizon [HUE!]. We need to solve (13.361) with the condition 

X(u,p) p- 1+A -J(u) asp^O (3.38) 

where A± are defined in (I3.28p . To perform the adiabatic expansion, let us decompose 
the field x{Pi u ) as 

X(p,u) = Xi(p,u) + Xs(p,u) (3.39) 

Where xi(.Pi u ) — J(u)p~ 1+A ~ and Xs(p,u) ~ p~ 1+A+ as p — > 0. 

For a constant J, xi(p,u) = Xi(p) is time independent. In this case there is a static 
solution Xs(p,u) = Xo(p) ; which is the equlibrium configuration. In the presence of a 
source which is slowly varying in units of the horizon radius, one can therefore expand 
the field Xs{p,u) in an adiabatic expansion of the form 

Xs(p,u) = Xo(p,J(u)) + e Xi(p,u) + ■■■ . (3.40) 

Here e ~ d u (recall that we are using r = 1 units) is an adiabaticity parameter which 
keeps track of the adiabatic expansion. If we scale u-fu/f, each u derivative is of order 
0(e). The idea then is to insert (13 .401) into the equations of motion and obtain equations 
for xii Xii ' ' ' order by order in e. To the lowest order one gets 

Vfxi = {[~d 2 p + V (p)} + f(r)(3xl + ^XiXo + ^)} X i = ~2d u d pX i - 2d u d pX o (3.41) 
The solution to this equation is 

Xx = I dp'G(p,p')d u ,d p ,( Xo + Xi){p')- (3-42) 



where G(p, p') is the Green's function of the operator T> p 1 ^ with the boundary conditions 
G(0,p) = and G(oo,p) is regular : 

g (1 W) = 3 g : ii(p'Mp), P<P' 

1 Up%(p), p>p', (3.43) 



where and ^2 are solutions of homogeneous part of eqn (I3.42p satisfying appropriate 
boundary condition at the horizon p = 00 and the boundary p = respectively, and 
W(^i,^2) is the Wronskian which is independent of p in this case. We have normalized 
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£i(p) and ^(p) i n such a fashion that £x = 1 at the horizon and £2 — > p~ 1+A ~ near the 
boundary. Regularity of the functions xi and Xo mean that d r xi,d r Xo are finite at the 
horizon. Since (r — 1) ~ e~ p , this implies that d p (xi + Xo) ~ exp(— p). This ensures that 
the integral in (13.421) is finite even though the Green's function approaches a constant in 
the region near the horizon p' — > 00. Furthermore, near the horizon £ 2 can be expressed 
as a linear combination of a regular and irregular solution, i.e ^(p — > 00) = ap + b. This 
implies that W(£i,f;2) = a - Thus Xi(w,p) is finite so long as a is finite. 

At the critical point, J becomes small. Then xo and xi hi t ne left hand side of 
(I3.4ip vanish, and the operator is identical to the operator acting on the linearized small 
fluctuations at m 2 = m 2 around the trivial solution xo — 0, i.e. the operator V p which 
appears on the right hand side of f 1 3 . 3 j) . We know that this operator has a zero mode 
which is regular at the horizon and vanishes as p~ 1+A ~ at the boundary p = 0. This 
means that at this point a = 0. Therefore, the first adiabatic correction diverges. For 
small J(u), the leading departure from the critical operator comes from the term which is 
proportional to Xo ~ J 2 ^ 3 - Thus we can use perturbation theory in J to estimate a oc J 2 / 3 . 
As argued before, xo ~ ( — J)K while xi ~ J- Hence the leading divergence in xi can be 
estimated as Xi{Pi u ) ~ J~ A ^J- Adiabaticity breaks down when 

Xl (p,u)~Xo^ J~ J 5/S (3.44) 

In particular for profiles of J(u) where J(u) ~ vu near the critical point at J = 0, 
adiabaticity breakdown occurs at 

u ~ v~ 2/5 (3.45) 



3.4 Quench in a Landau Ginsburg Model 



The scaling behavior found above is identical to that in a Landau- Ginsburg dynamics 
with dynamical critical exponent z = 2. The dynamics of an spatially homogeneous order 
parameter <p is given by 

% + m 2 ip + cp 3 + J(t) = (3.46) 
at 

The equilibrium critical point is at m = J = 0. For m 2 = the equilibrium value of the 
order parameter is 

tpM = [-J] 1 ' 3 (3.47) 
As usual an adiabatic expansion is of the form 

<p(t) = <p {J(t)) + ei Pi (*) + ••• (3-48) 



16 



and to lowest order 



<Pi 



1 -d(p 

r J ~ 



2ipl dJ 

and adiabaticity breaks down when ipx ~ ip which becomes the condition 



J 



dJ 



JV3 ^ j „ J5/3 



(3.49) 



(3.50) 



exactly as in our system. 

When adiabaticity breaks down our system enters a scaling region. Suppose the func- 
tion J(t) behaves linearly with time in the critical region. Then it is straightforward to 
see from ( I3.46P that in this region the solution is of the form 

<p(t,v) = v 1/5 V (tv 2/5 ,l) (3.51) 

This means, in particular, that the time at which the order parameter hits zero scales as 
i/~ 2 / 5 while the value of the order parameter at t — scales as v 1//5 . A numerical solution 
of the equation (I3.46P with adiabatic initial conditions is consistent with this scaling, as 
shown in Figure ([8]) 




Figure 8: The scaled order parameter as a function of scaled time for a J(t) = tanh(t> t) 
at m 2 = with v = 10~ a5 , 10 _1 , 10~ L5 , 10" 2 (from the bottom on the left). The adiabatic 
solution (dashed) is also shown as a comparison. 

Note that the order parameter hits zero later than the location of the equilibrium 
critical point. This is a manifestation of the phenomenon of raising the critical tempera- 
ture when the temperature is time dependent (35] which has been holographically realized 
in [36] . 

We will now argue that the behavior of our holographic system in the critical region 
is fairly well described by such a LG dynamics 
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3.5 Small v dynamics in the Holographic Model 



Consider the dynamics of the bulk field in the critical region in the presence of a linear 
quench J(u) = vu for small v. To do this, first substitute (I3.39P in the equation ( I3.36P 
and rescale 

Xs -> v^Xs, u v~%u (3.52) 
The equation f!3.36j) then becomes 

[-0? + V (p)}xs + vl[f(r)(xs) 3 + u[-d 2 p + V (p)]xi + 2d a d p x s } + ■ ■ • = (3.53) 

The ellipsis denote terms which contains higher powers of v. 

Let us expand the sub-leading part of the scalar field in terms of eigenfunctions of the 
operator V p (defined in equation (I3.30P at the critical point, 



Xs(p,u) = J a k (u)xk(p)dk (3.54) 

The Xk satisfy 

V c p Xk = [~d 2 p + V c (p)] Xk = k\ k (3.55) 

where V C denotes the potential in (13.31!) at m 2 = m 2 . The eigenfunctions Xfc(p) are delta 
function normalized and obey the condition 

Um^oip^-Xkip)} = (3.56) 
In terms of the eigen-coefficients a k {u) the equation f!3.53|) becomes, 

k 2 dk + V s ^Mj7fc + J bkk'duO-k'dk' + J ayayiaynCk^ 1 \k" ,k"'dk ' dk' 'dk'j + • • ■ = (3.57) 



where 



Jk = j Xk(p)[-d 2 p + V (p)]xidp 
hk> = J dpxkdpx'k 
C k y,k",k"' = / dpXkXk'Xk"Xk"'f(r). (3.58) 



The equation (I3.57P suggests that there is a solution in a perturbation expansion of powers 
of v i, 

a k {u) = 5{k)^ {u) + vlf} k {u) + ■■■ , (3.59) 
18 



If this expansion makes sense, the dominant behavior of the solution for Xs is given by 
the zero mode £ - 

Substituting ( 13.591) in the equation ( 13.571) we get to the lowest order in the small v 
expansion, 

uJo + b 00 ^o(u) + C ooo£o(m) 3 = 

Vk{u) = -^j (uj k - b ko-^io(u) - Cfcooofo(M) 3 ^ (3.60) 

The first equation in f l3.60f) determines the time dependence of the zero mode £o ; while 
the second equation determines the leading correction from nonzero modes in terms of 
the solution for £ - 

It is useful to rewrite the second equation by subtracting the first from it, 

Vk(u) = --^ (u(Jk - Jo) - (b k0 - & oo)^t6d(m) - (Cfeooo - Coooo)fo(w) 3 ^ (3.61) 

The expansion in powers of v 2 ^ would be valid if r)p.(u) remains finite. However, k is a 
continuous parameter starting from zero. This means that there is a potential divergence 
in the k — > limit. 

Indeed, as will be argued in the next subsection, for generic potential Vq the numerator 
on the right hand side of (13.611) behaves as k for small k, so that % indeed diverges at 
k = 0. However exactly at the critical point, the small k behavior changes to k 2 so that 
rjk remains finite and the expansion in v 2//5 remains valid. 



3.6 Validity of the small v expansion 



To examine this issue we need to consider the eigenvalue problem 

[~d 2 p + V Q (p)]Xk = k 2 X k (3.62) 

As discussed above the potential Vo(p) — > —e~ p as p — > oo. This potential is shown in 
Figure ([7]). 

The basic features of the eigenfunctions can in fact be gleaned from a simpler problem 
in which we replace the potential by the following potential which has the same qualitative 
features. 

= oo p = 

U{p) = -U 0<p<l 

= 1 < p < oo (3.63) 
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This problem is of course solvable. The eigenfunctions of the Schrodinger operator with 
eigenvalue k 2 > are 



MP) = ^£sm(v^TT^p) 

\/7T 



sin(kp + 8(p)) 



< p < 1 

< p < oo 



(3.64) 



where the constants A(k) and 8(k) are determined by matching at p = 1, 
A(k) = ' 



k 2 sm 2 (VW+Uo) + Vk 2 + U cos 2 Wk 2 + U ) 



9(k) 



tan 



A; tan (y/k 2 + U Q ) 
Vk 2 + U ~ 



-k. 



(3.65) 



The solution for k = is 



B 



71 



sin(v / t/o P) 



< p < 1 



ap + b 



< p < oo 



(3.66) 



The matching conditions at p = 1 now yield 

B 



—=sm(\/Uop) = a + b 

\/7T 



cos(v / t/o p) 



5^ 



(3.67) 



For any a ^ the solution blows up at p = oo. Thus regular solutions require a = 0. 
However the second equation in ( 13.67!) then imply that 



(3.68) 



These are the zero modes. As we increase the depth of the potential, the first zero mode 
appears at \/Uo = tt/2. In the context of our model this is the potential where we have a 
critical point. 

The small k behavior of A(k) and 8(k) can be read off from the expressions (13.651) . 
For a generic Uq these are 

k 



A(k) 



/ Uq COS v^O 



0(F 



W „ fc[ W&o _ l] + 0(A: 3 ) 
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(3.69) 



whereas for critical potentials we have 

A(k) ~ l- — + 0{k 



9(k) ~ (3.70) 

Thus the small-fc behavior of the eigenfunctions are drastically different for the critical 
potentials. This has important implications for the coefficients like (Jk — Jo)-, (b^o — 600) 
and (Cfcooo — Coooo) in fl3.61[) . Consider for example the quantity Jk- This is an integral 
of the form 

/■oo 

dp J{p) Xk(p) (3.71) 



where J(p) is a smooth function (which is [—d 2 p + Vo(p)]xi)- If we replace the true eigen- 
functions by those of our simplified problem, we get 

/1 i* 00 

dp sin(V 'k 2 + U p) J (p) + J dpsm(kp + 6(k)) (3.72) 

Using (I3.69P and ( I3.70P we therefore see that 

Jk-Jo-k k ->• (3.73) 

for generic potentials, whereas 

Jk-Jo-k 2 k-+0 (3.74) 

for critical potentials. It is straightforward to see that the behavior of the other coefficients 
(ho ~ b 00 ) and (C k0 oo - C 00 oo) are similar. 

The small-fc behavior of the eigenfunctions for the potential which is relevant for us. 
V Q c (p) is quite similar. This has been discussed in the Appendix C of [TT| . Going back 
to (I3.6ip we therefore see that the small v expansion is generically not valid since the 
corrections diverge at small k. However for the critical potential, fj^ remain finite as 
k — > and the expansion in powers of v 2//5 makes sense. 

3.7 Scaling relation 

The validity of this expansion means that the solution for Xs is of the form 

1 ~ 2 3 / 2 

Xs(p,u) w vs£o(vsu)xo(p) +v f ' / fjk(v*u)xk(p)dk-\ (3.75) 
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Thus the dynamics for small v is dominated by the zero mode in the critical region. The 
equation for £ is, however, exactly the same as the equation for the order parameter <p in 
the LG model in the previous subsection, after going back to the original variables prior 
to the rescaling in (I3.52p . Thus in this region the system behaves as one with dynamical 
critical exponent z = 2. Therefore we conclude that to leading order in small v , we have 

< O > (it, v) ~ Lim^o[p 1_A+ X s (p, u)] ~ v 1 ' 5 <0> (tv 2 '\ 1) (3.76) 

and the time scale behaves as v~ 2 / 5 . 

The fact that the dynamics in the critical region is governed by an equation with a 
first order time derivative is made manifest in our treatment using Eddington-Finkelstein 
coordinates. The whole analysis can be of course performed in principle in the (t, r) 
coordinates. However, we suspect that in this case one has to exercise extreme care, just 
as one had to extract out a leading horizon behavior in the linearized problem of fields in 
a black brane background [37] . 

Note that beyond the critical region, J(u) departs from the form J(u) ~ vu, and 
the expansion in powers of v 2 ^ is no longer valid. Now all the modes are important, 
and integrating out higher modes can give rise to higher time derivatives in the effective 
equation for < O >. 

3.8 Mass quench 

An analysis similar to the above can be carried out when the quenching is performed 
by making the mass parameter of the bulk theory a function of the retarded time u, 
keeping J(u) = 0. This does not have a direct interpretation in the boundary field theory. 
However, as explained in [TjJ], the field can acquire a mass because of a coupling to some 
other field <p' . A time dependent boundary value of the field 0' can then lead to a time 
dependent mass. When m 2 (u) ~ m 2 + vu near the critical point, we now get a behavior 
< O >~ v x l\ 

3.9 Numerical Results 

We have performed some preliminary numerical work for the case of a mass quench. Our 
results clearly display the propagation of disturbances towards the horizon along a light 
cone , and a decay of the order parameter in a fashion similar to the z = 2 LG dynamics, 
which is shown in Figure (jHJ). However our results are not accurate enough to verify the 
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scaling behavior found above. The late time decay should be governed by the quasinormal 
modes [3l|39]. 



<0>(t) 
0.010 f- 




-2 2 4 6 

Figure 9: The order parameter < O > (t) in boundary theory. 

4 Outlook 

We have shown that holographic methods are useful in providing insight into various 
questions related to quantum quench in strongly coupled field theories which have a 
gravity dual. Interestingly, this is possible in the probe approximation, which is typically 
much easier to study. 

So far we have been able to study in some detail quench dynamics near holographic 
critical points which are described by mean field exponents. Not surprisingly, we found 
scaling behavior characteristic of Landau-Ginsburg models with z — 2. This value of the 
dynamical critical exponent is consistent with the results of [38J. 

It is important to perform our quench analysis for zero temperature, where the equi- 
librium transition is of the BKT type [19], similar to brane models of zero temperature 
chiral symmetry breaking transition in [32J. This case is rather subtle, but can be studied 
using similar methods. We expect that in this case we will have a z = 1 theory dominating 
the critical region. This would provide results for quench dynamics which are not easily 
obtainable by other methods. Finally, a more extensive numerical investigation should 
throw light on the question of thermalization at late times, after the system has crossed 
the critical region. 
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